Analytical solution to position dependent mass Schrodinger equation 
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Using a recently developed technique to solve Schrodinger equation for constant mass, we stud- 
ied the regime in which mass varies with position i.e position dependent mass Schrodinger equa- 
tion(PDMSE). We obtained an analytical solution for the PDMSE and applied our approach to 
study a position dependent mass m(x) particle scattered by a potential V(x). We also studied the 
structural analogy between PDMSE and two-level atomic system interacting with a classical field. 
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Schrodinger equation with position dependent mass is 
one of the areas of research which has gained great at- 
tention in the past decades. Position-dependent mass 
Schrodinger equation(PDMSE) has been applied to sev- 
eral physical systems. For example PDMSE is applied 
in electronic properties of semiconductors pQ, quantum 
dots and quantum wells [21 E], semiconductors hetero- 
structures @], supper-lattice band structures [5], He- 
Clusters [5] quantum liquids [7], the dependence of en- 
ergy gap on magnetic field in semiconductor nano-scale 
rings [5], the solid state problem with Dirac equation [9 
etc. One of the motivation behind investigating these sys- 
tem with position dependent mass is, how do these mass 
variation effects the dynamics of the quantum-mechanical 
system. In such systems the energy is described by a 
Hamiltonian which contains the kinetic energy and the 
potential energy opeators Jt? = 2? + Y, where special 
care is taken for the kinetic term. 

O. von Roos jTU] was the first to suggest the follow- 
ing generalized form of the kinetic energy operator for 
position-dependent mass model 

3" = -(m r, pm e pm p + m p pm e pm n ) ) (1) 

where m = m(r) is the position-dependent mass. The 
constants 77, e and p, which are also knows as the von 
Roos ambiguity parameters can be assumed to be arbi- 
trary but they obey the constraint equation 77 + e + p = 
— 1. Before von Roos several forms of operator ST has 
been used to solve this problem [TTH14] . Different ap- 
proaches have been used to find analytical solution to 
PDMSE like point-canonical transformation |15| . Green's 
function [16] , Heun equation [TTl [18] , Group-Theoretical 
method Q35], potential algebra [2U], Lie-algebraic [3T] and 
supersymmetric [S3] approach etc. 

In this Brief Report, we obtained an analytical solu- 
tion for the position dependent mass Schrodinger equa- 
tion with general mass variation m(x) by transforming 
the PDMSE to Riccati equation. Analytical solution 
for Schrodinger equation with constant mass, beyond 
adiabatic approximation, has been recently investigated 
extensively [53] where we applied our approach to 1-D 



[2"3"] and 3-D [M] scattering problem and showed that 
our method gives better accuracy than the well-known 
JWKB method. Here in this article we have extended 
that approach to variable mass regime and obtained very 
accurate results for PDMSE. The main result of this pa- 
per is the analytical solution given by Eq.(19|. To illus- 



trate how well our method works we considered a posi- 
tion dependent mass m(x) particle scattered by potential 
V(x) and obtained the wave function for ID case (which 
could easily be extended to 3D). Plots of the numerical 
simulation and the analytical solution are nearly identical 
[see figs 1,2,3]. We also briefly discuss structural anal- 
ogy between PDMSE and two-level atom(TLA) driven 
by resonant classical field thus bridging between quan- 
tum optics and condensed matter Physics. 



I. MODEL: ANALYTICAL SOLUTION 

Let us consider the symmetric ordering form of the ki- 
netic energy operator as given by Eq.([T]) i.e (77 = 0, e = 
— 1, p = 0). The Hamiltonian of the system is written as 



1 

P— P 

m 



r(r). 



(2) 



In 1-D, this Hamiltonian gives the following form of the 
Schrodinger equation with position dependent mass, 



y dx 



1 (Nf(x) 



m(x) dx 



+ V(x)*(jc) =£V{x). (3) 



Here the mass function m{x) — m Q g(x). Let us non di- 
mensionlize Eq. ^ using the following scaling parameters 



E 



S/£ , V(x) = V(x)/£ , z = (V2m fo/fi) x. 



we get, 



d 

dz 



1 d&(z) 

q{z) dz 



(4) 



+ V(x)^(z)=E^(z). (5) 



2 



Let us recast Eq.([5]) in a more desirable form we need, Let us consider ei(^o) ~ which gives B = hence 



f-2 ~ + e(*)[E - V(z)]\ *(*) = 0. (6) 

dz z g(z) dz J 

Let us define two functions and £(z) as, 
2g{z) 

az) = g(z)[E-V(z)). (7b) 
Substituting Eq.([7) in Eq.Q we get, 

To find an analytical solution for the Eq.([8]) we make a 
formal substitution of 



9(z) 



(9) 



which reduces the Schrodinger equation to non-linear 
Riccati equation of the form 



/'(z) + /(z) 2 -2C(z)/(z)+eW=0. 



(10) 



Let us work in the adiabatic regime of the problem i.c to 
the zeroth order and neglect the contribution of f'(z) in 
Eq.(pJ). Solving Eq.dlOjl for f{z) = we get, 



(11) 



The zeroth order general solution of the Schrodinger 
equation with position dependent mass Eq.j8j is 



*(z) = A ie ^ z)+l6(z) + A 2 e^ z) - ie{z \ 



where, 



ip(z) 



({z)dz, 



(12) 

(13a) 
(13b) 

To go beyond the zeroth order approximation assume, 
h(z) = f (z) + e 1 (z). (14) 

Substituting Ec . 



0(z) = / v/£(2) - ((z)*dz. 



cx e (z), Eq.(10l gives 



(14 1 in Eq.(10) and neglecting the term 
gives 

f Q {z) + e'^z) + 2/ (z) £l (z) - 2C(z)ei(z) = 0. (15) 



The general solution to Eq( 15 ) is given as 

&(z)e-^ 



where, 



dz + B } e' 



(16) 



u{z)=2 / [-f (z") + ((z")]dz" 

J z 



7o(2)e 



-u(z) 



dz> e 



(18) 



The first order general solution of the Schrodinger equa- 
tion with position dependent mass is 



= Axexp[(p(z) +i6(z) - <p + (z)] 
+A 2 exp [tp(z) - i0(z) - (f>-(z)] , 



(19) 



where f± (z), ip(z), 0{z) is given by Eq.( 11 ), Eq.( 13 i) and 



Eq.(13}) respectively and <fi±(z) is defined as 

<t>±{*) = ri r \f ± (z)e^ z >-^ dz 



dz. (20) 



Using an iterative procedure our approach can be easily 
extended to include the next order corrections. Indeed, 
assuming /2(z) = fo(z)+ei{z)+€2(z) we can easily obtain 
the equation for 62(2) as 

e' 2 (z) + 2 [f (z) + aiz) - C(z)] e 2 (z) + e\{z) = 0. (21) 

In general, for f n (z) — fo{z) + X)j=i e j( z )i we can obtain 

n-1 



e' n (z) + 2 



Mz) + ^ £j -C(z) 



e n (z) + e 2 n _ 1 (z)=0. 

(22) 

All these equations can be solved exactly, thus we have 
obtained an analytical approximate solution for PDMSE. 



A. Application to scattering problem 

Let us now compare our approximate analytical solution 
with the exact solution for some physical problem. As an 




FIG. 1: (Color Online) Gaussian potential V(z) = 
exp [— /3 2 z 2 ] and the hyperbolic mass variation m(z) = (mi + 
m 2 )/2 + [(mi - m 2 )/2] tanh(az) [see inset]. Plot of |*(a)| 2 
versus z for the energy E=2.5. Dashed line represents the 
exact solution and solid line is the solution from Eq.( 19 1. For 



(17) 

numerical simulation m 2 = 2, mi = 1, a = 0.1 and 7 = 0.5 
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FIG. 2: (Color Online) Gaussian potential V(z) = 
cxp [— /3 2 z 2 ] and the Gaussian-mass variation m(z) = 1 + 
exp [— 5 2 z 2 ] [see inset]. Plot shows ^(z)! 2 versus z for the 
energy E=1.75. Dashed line represents the exact solution and 
solid line is the solution from Eq.(19l. For numerical simula- 
tion /3 = 1/(2^5), 5 = 0.1 



FIG. 3: (Color Online). Gaussian potential V(z) = 
— sech(7z) and the mass variation m(z) — (mi + m2)/2 + 
[(mi — m,2)/2] tanh(az) [see inset]. |'I'(2)| 2 versus z for the 
energy E=0.1. Dashed line represents the exact solution and 
solid line is the solution from Eq.(19l. For numerical simula- 
tion wi2 = 2, mi = 1, a = 0.1 and 7 = 0.5. 



example, we study the scattering of a particle by a 1-D 
potential. A particle is propagating from left to right, 
during the propagation the effective mass of the particle 
varied due to its interaction with an ensemble of particles. 
For numerical simulations, we have used the scaling pa- 
rameters given by Eq.Q and the corresponding PDMSE 
Eq.(j8|. From elementary quantum mechanics we know 
that due to the interaction of the particle with the po- 
tential, for z — > +00 the wave function for the particle is 
a plane wave while for z — > — 00 it is the sum of the in- 
cident and the reflected plane waves. Thus we can write 



the mass function from rri\ to mi [see Figs. (1,3) inset] 
can be controlled by the parameter a or S as appropriate. 
Last example which we considered for the potential func- 
tion is V(z) = — sech(7z) and the plot of the probability 
density (both numerical and analytical) is shown in Fig 
3. Numerical simulation and analytical solution shown in 
Figs. (1,2,3) are nearly identical which suggests that our 
method allows us to accurately obtain the wave function 
of a position dependent mass m{x) particle scattered by 
a potential V(x). Our approach can be easily extended 
to find bound states also. 



*(z) -> A$i(z) + BV r (z) 



(23) 



Here the incident and the reflected plane wave have the 
form 



*,•(*) 
* r (z) 



exp 
exp 



WEz 
-i4~Ez 



(24a) 
(24b) 



As a first example we will consider a Gaussian potential 
and the hyperbolic mass function of the for 

V{z) = exp [~fi 2 z 2 ] , (25a) 

\m 1 +m 2 ] [mi-m 2 ] , 

m(z) — + tanh(az). (25b) 

The mass m(z) of the particle changes from mi — > 
Fig. (1) shows the plot of the the probability density 
|\l/(z)| 2 against z. The dashed line is the exact solution 
of the PDMSE Eq.Q while solid line is the analytical 
solution Eq. (Il9 1 . 



Keeping the same form for the potential let us takes 
the mass function as m(z) = 1 + exp [— 5 2 z 2 ~\ . Here Fig. 
(2) shows the comparison plot for the numerical simula- 
tion and our analytical result. The smooth variation of 



B. Position Dependent Mass Schrodinger Equation 
and Two-Level System 

Let us consider a TLA with |a) and \b) representing the 
upper and the lower level states [see Fig. 4(a)]. The 
equation of motion for the probability amplitudes C a and 
Cb can be written as [5S] 

C a (t) = in(t)cos{vt)e iut C b (t), (26a) 
C b {t) = in(t)cos(vt)e- iut C a (t). (26b) 

where Huj is the energy difference between two levels, f2(f ) 
is the Rabi Frequency and considered real here. In invok- 
ing rotating-wave approximation we let cos(ft)e ±lult — > 
e ±lA */2 where A = cj — v 1 is detuning from resonance. 

(27a) 
(27b) 



C a = i p(t)/2] e lAt C b , 
C b - i p)(t)/2] e~ lAt C a 



Coupled first order equations Eq.(27) can be transformed 
to Riccati equation for / = C a (t)/C b (t) [26l[27] or to the 
second order differential equation [2H1 HI] 



C a (t) ~ 



iA + 



rut) 



c a (t) 



Sl 2 (t) 



C a {t)=0. (28) 
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Solution to Eq. ( 32 ) is given as 



FIG. 4: (Color Online) (a) Two-level atomic system, atomic 
transition frequency uj — ui a — U)i, detuning A = u> — v and 
Rabi frequency Q(t). (b) Few cycle sine (dashed line) and 
cosine (solid line) pulse with gaussian envelope. 



In the case of resonance uj = v these equations Eq.(28) 
transforms to 



C a (t) 



C a {t) 



n 2 (t) 



C a {t) = 0. (29) 



The exact solution for the transformed equation is given 
as 



C a (t) =A ieX p 
A 2 exp 



dt'Sl{t')/2 



dt'Q(t')/2 



(30) 



Here A\ and A 2 are constant which can be obtained 
using the initial condition of the problem. To see the 
connection between PDMSE and two-level system (TLS) 
let us consider that the mass and potential function is 
related as 



g(z) + V{z) = E 
Subsequently from Eq.(|6]) we get 
f d 2 g'(z) d 



dz 2 q[z) dz 



(31) 



(32) 



= B ie ^ (2) +B 2 e- ltp{z \ 
•f{z) = \ g{z)dz. 



(33a) 
(33b) 



If we look Eq.(32| with Eq.(29) we see that they share 



the same structure. In one case it is the Rabi frequency 
which drives the TLA, while the variable mass g(z) 
plays a similar role for PDMSE. Thus the equivalence 
can be summarized as 



ZHt, 

o C a (t), 

e {z) <-> n(t)/2. 



(34a) 
(34b) 
(34c) 



Here we have assume that Q(t) remains positive for all 
t. This structural analogy between the TLS and the 
PDMSE motivates us to find similar connections between 
two quantum-mechanical systems in two different areas 
of Physics. 

To conclude, in this article we have studied the dy- 
namics of position-dependent mass Schrodinger equation 
using a recently developed technique (for constant mass 
Schrodinger equation) to find analytical solution. The 
approximate solution is in excellent agreement with the 
numerical simulations [see Figs. 1,2,3]. We also investi- 
gated the structural analogy between coherent excitation 
of a TLA with a classical field in RWA and PDMSE. 
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